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Abstract 

What is the higher-dimensional analog of a permutation? If we 
think of a permutation as given by a permutation matrix, then the 
following definition suggests itself: A d-dimensional permutation of 
order n is an n x n x . . . n = [n] d+1 array of zeros and ones in which 
every line contains a unique 1 entry. A line here is a set of entries of 
the form {(x±, . . . , Xi_i,y, x^+i, . . . , Xd+i)\n > y > 1} for some index 
d + 1 > i > 1 and some choice of Xj 6 [n] for all j ^ i. It is easy to 
observe that a one-dimensional permutation is simply a permutation 
matrix and that a two-dimensional permutation is synonymous with 
an order-n Latin square. We seek an estimate for the number of d- 
dimensional permutations. Our main result is the following upper 
bound on their number 



n - " d 



+ 

We tend to believe that this is actually the correct number, but the 
problem of proving the complementary lower bound remains open. 
Our main tool is an adaptation of Bregman's [1] proof of the Mine 
conjecture on permanents. More concretely, our approach is very 
close in spirit to Schrijver's [11] and Radhakrishnan's [10] proofs of 
Bregman's theorem. 
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1 Introduction 



The permanent of an n x n matrix A = (a^) is defined by 

n 

Per(A)=Y,U a ^ 

aeS n i=l 

Permanents have attracted a lot of attention [9]. They play an important 
role in combinatorics. Thus if A is a 0-1 matrix, then Per(A) counts perfect 
matchings in the bipartite graph whose adjacency matrix is A. They are 
also of great interest from the computational perspective. It is #P-hard 
to calculate the permanent of a given 0-1 matrix [12], and following a long 
line of research, an approximation scheme was found [6] for the permanents 
of nonnegative matrices. Bounds on permanents have also been studied at 
great depth. Van der Waerden conjectured that Per(A) > for every 
n x n doubly stochastic matrix A, and this was established more than fifty 
years later by Falikman and by Egorychev [4, 3]. More recently, Gurvitz [5] 
discovered a new conceptual proof for this conjecture (see [8] for a very 
readable presentation). What is more relevant for us here are upper bounds 
on permanents. These are the subject of Mine's conjecture which was proved 
by Bregman. 

Theorem 1.1. If A is an nx n 0-1 matrix with ones in the i-th row, then 

n 

Per (A) < H(^!) 1/ri - 

i=l 

In the next section we review Radhakrishnan's proof, which uses the en- 
tropy method. Our plan is to imitate this proof for a <i-dimensional analogue 
of the permanent. To this end we need the notion of d- dimensional permu- 
tations. 

Definition 1.2. 1. Let A be an [n] d array. A line of A is vector of the 
form 

(A(ix, ...,i i _ 1 ,t,i j+ i J ...,z d ))t =1 , 
where 1 < j < d and i\ , . . ., ij-i, , . . . , id G [n] . 
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2. A d- dimensional permutation of order n is an [n] d+1 array P of zeros 
and ones such that every line of P contains a single one and n — 1 
zeros. Denote the set of all d-dimensional permutations of order n by 



For example, a two dimensional array is a matrix. It has two kinds of 
lines, usually called rows and columns. Thus a 1-permutation is an n x n 0-1 
matrix with a single one in each row and a single one in each column, namely 
a permutation matrix. A 2-permutation is identical to a Latin square and 
S2 n is the same as the set C n , of order-ro Latin squares. We now explain the 
correspondence between the two sets. If X is a 2-permutation of order n, 
then we associate with it a Latin square L, where L(i,j) as the (unique) index 
of a 1 entry in the line A(i,j,*). For more on the subject of Latin squares, 
see [13]. The same definition yields a one-to-one correspondence between 
3-dimensional permutations and Latin cubes. In general, d- dimensional per- 
mutations are synonymous with <i-dimensional Latin hypercubes. For more 
on d- dimensional Latin hypercubes, see [14]. To summarize, the following is 
an equivalent definition of a <i-dimensional permutation. It is an [n] d array 
with entries from [n] in which every line contains each i £ [n] exactly once. 
We interchange freely between these two definitions according to context. 

Our main concern here is to estimate |*Srf iri |, the number of d-dimensional 
permutations of order n. By Stirling's formula 



As we saw, |S2. n | is the number of order n Latin squares. The best known 
estimate [13] is 



This relation is proved using bounds on permanents. Bregman's theorem for 
the upper bound, and the Falikman-Egorychev theorem for the lower bound. 
This suggests 

Conjecture 1.3. 







In this paper we prove the upper bound 
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Theorem 1.4. 

l^»|<((l + o(l)) ? ) . 

As mentioned, our method of proof is an adaptation of [10]. We first need 

Definition 1.5. 1. An [n] d+1 0-1 array M\ is said to support an array 
M 2 if 

M 2 (i u ...,i d+1 ) = 1 M x {i u ...,i d+1 ) = 1. 

2. The d-permanent of a [n] d+1 0-1 array A is 

Per d (A) = The number of d- dimensional permutations supported by A. 

Note that in the one- dimensional case, this is indeed the usual definition 
of Per(A). It is not hard to see that for d = 1 following theorem coincides 
with Bregman's theorem. 

Theorem 1.6. Define the function f : N>o x N — > R recursively by: 

• /(0,r) = log(r), where the logarithm is in base e. 

• Mr) = JEU /(<*-!,*). 

Let A be an [n] d+1 0-1 array with r^...^ ones in the line A(ii, id, *)■ Then 

Per d (A) < Y[ e /(^n,-- d ). 

h,...,id 

We will derive below fairly tight bounds on the function / that appears 
in theorem 1.6. It is then an easy matter to prove theorem 1.4 by applying 
theorem 1.6 to the all-ones array. 

What about proving a matching lower bound on Sd, n (and thus proving 
conjecture 1.3)? In order to follow the footsteps of [13], we would need a lower 
bound on Per d A, namely, a higher-dimensional analog of the van der Waer- 
den conjecture. The entries of a multi-stochastic array are nonnegative reals 
and the sum of entries along every line is 1. This is the higher-dimensional 
counterpart of a doubly-stochastic matrix. It should be clear how to extend 
the notion of Per d (A) to real-valued arrays. In this approach we would need 
a lower bound on Per d (A) that holds for every multi-stochastic array A. 
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However, this attempt (or at least its most simplistic version) is bound to 
fail. An easy consequence of Hall's theorem says that a 0-1 matrix in which 
every line or column contains the same (positive) number of 1-entries, has 
a positive permanent. (We still do not know exactly how small such a per- 
manent can be, see [8] for more on this). However, the higher dimensional 
analog of this is simply incorrect. There exist multi-stochastic arrays whose 
(/-permanent vanishes, as can easily be deduced e.g., from [7]. 

We can, however, derive a lower bound of \Sd, n \ > exp(Q(n d )) for even n. 
Consider the following construction: Let n be an even integer, and let P be 
a (/-dimensional permutation of order [|] . It is easy to see that such a P 
exists. Simply set 

P(i 1 ,...,i d ) = (i 1 + ... + i d ) mod-. 

Now we construct a (/-dimensional permutation Q of order [n] d by replacing 
each element of P with a [2] d block. If P(i\, id) = j, then the correspond- 
ing block contains the values j and j + | . It is easy to see that there are 
exactly two ways to arrange these values in each block, and that Q is indeed 
a (/-dimensional permutation of order [n] d . There are (^) blocks, and so the 

number of possible Q's is 2^) . For a constant d this is exp(f2(n d )). 

In section 2 we present Radhakrishnan's proof of the Bregman bound. 
In section 3 we prove theorem 1.6. In section 4 we use this bound to prove 
theorem 1.4. 

2 Radhakrishnan's proof of Bregman's theo- 
rem 

2.1 Entropy - Some basics 

We review the basic material concerning entropy that is used here and refer 
the reader to [2] for further information on the topic. 

Definition 2.1. The entropy of a discrete random variable X is given by 
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For random variables X and Y , the conditional entropy of X given Y is 

H{X\Y) = E[H(X\Y = y)] = ^Pr(F = y)H{X\Y = y). 

y 

In this paper we will always consider the base e entropy of X which simply 
means that the logarithm is in base e. 

Theorem 2.2. 1. If X is a discrete random variable, then 

H{X) < log \range(X)\, 
with equality iff' X has a uniform distribution. 

2. If Xi, X n is a sequence of random variables, then 

n 

H(X 1 , X n ) = 2^ H(Xi\Xi, Xi-i). 
i=i 

3. The inequality 

H{X\Y) < H(X\f(Y)) 

holds for every two discrete random variables X and Y and every real 
function /(•). 

The following is a general approach using entropy that is useful for a 
variety of approximate counting problems. Suppose that we need to estimate 
the cardinality of some set S. If X is a random variable which takes values 
in S under the uniform distribution on S, then H(X) = logdS"!). So, a good 
estimate on H(X) yields bounds on \S\. 

This approach is the main idea of both Radhakrishnan's proof and our 
work. 

2.2 Radhakrishnan's proof 

Let A be an n x n 0-1 matrix with r,j ones in the i-th row. Our aim is to 
prove the upper bound 



Per(A) < ] 



i=l 
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Let M. be the set of permutation matrices supported by A, and let X be 
a uniformly sampled random element of M.. Our plan is to evaluate H(X) 
using the chain rule and estimate using the fact (theorem 2.2) that 
H(X) = \og(\M\). 

Let Xi be the unique index j such that X(i, j) — 1. We consider a process 
where we scan the rows of X in sequence and estimate H(X) = H(Xi, X n ) 
using the chain rule in the corresponding order. To carry out this plan, we 
need to bound the contribution of the term involving Xi conditioned on the 
previously observed rows. That is, we write 

n 

H(X) = ^H(X i \X 1 ,...,X i - 1 ). 

Let Ri be the set of indices of the 1-entries in A's f-th row. That is, 

Ri = {j:A(i,j) = l}. 

Let 

Zi = {j G Ri : X^ = j for some i' < %}. 

Note that Xi e because X is supported by A. In addition, given that 
we have already exposed the values X^ for z' < z, it is impossible for Xi to 
take any value j G Zi, or else the column X(*,j) contains more than a single 
1-entry. Therefore, given the variables that precede it, Xi must take a value 
in Ri \ Zi. The cardinality iVj = \Ri \ Zi\ is a function of X±, and 
so by theorem 2.2, 

n 

H{X) = Y i H(X i \X 1) ...,X i . 1 ) 

n 

= ^ ^ Pr(Xi = xi, ...,Xi-x = Xi-i)H(Xi\Xi = x x , = a;»_i) 

i=l xi,...,Xi—i 

n 

<Y, Yl P<Xi=Xi,..;X^ 1 = x l . 1 )\og{N i ) 

1 = 1 Xl,...,Xi—l 

n n 

= Yl ^ v [log(iVi)] = $]E X [logC^)]. 

i=l i=l 
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It is not clear how we should proceed from here, for how can we bound 
log(iVj) for a general matrix? Moreover, different orderings of the rows will 
give different bounds. We use this fact to our advantage and consider the 
expectation of this bound over all possible orderings. Associated with a 
permutation a £ § n is an ordering of the rows where Xj is revealed before 
Xi if a(j) < cr(z). We redefine Zj and iVj to take the ordering a into account. 
Let 

Zj(cr) = {j £ Ri : Xi> = j for some a(i') < cr(i)}. 

Ni(a) = \Ri^Zi(a)\. 

Then A^(cr) is the number of available values for Xi, given all the variables 
Xj for j such that a(j) < cr(i). As before, using the chain rule we obtain the 
inequality 

n n 

H{X) = Y^HiX^X, : a(j) < a({)) < [log(iV 4 (a))]. 

i=l i=l 

The inequality remains true if we take the expected value of both sides 
when a is a random permutation sampled from the uniform distribution on 

n n 

h(x) <J2 E ° tMW)]] = J2 Ex PosC^W)]]- 

1=1 1=1 
Thus, the bound we get on H(X) depends on the distribution of the random 
variable Ni(a). The final observation that we need is that the distribution 
of Ni(cr) is very simple and that it does not depend on X. Consequently we 
can eliminate the step of taking expectation with respect to the choice of X. 
Let us fix a specific X. 

Let Wi denote the set of r» — 1 row indices j ^ i for which Xj £ Ri. 
Note that iVj is equal to minus the number of indices in Wi that precede 
i in the random ordering a. Since a was chosen uniformly, this number 
is distributed uniformly in {0,...,r,; — 1}. Thus, iVj is uniform on the set 
{l,...,rj}. Therefore 

E a [logger))] = V - log(fc) = - login!). 



Hence 



71 r 1 -i n 1 

//(X)<^E X -log(rj) =^-log(r l !) 
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which implies the Bregman bound. 



3 The d-dimensional case 
3.1 An informal discussion 

The core of the above-described proof of the Bregman bound can be viewed 
as follows. Let us pick first a 1-permutation X that is contained in the matrix 
A and consider the set Ri of the r, 1-entries in A's i-th row. There are exactly 
Ti indices j for which Xj G Ri. The random ordering of the rows determines 
which of these will precede the i-th row (or will cast its shadow on the i-th 
row). The random number Ui of rows that cast a shadow on the i-th row is 
uniformly distributed in the range {0, . . . , Ti — 1}. The contribution of this 
row to the upper bound on H(X) is E CT [logiVj], where iVj = — Ui is the 
number of 1-entries in the i-th row that are still unshaded. The expectation 
of log Ni is exactly j- YJ'-=i lo Si = v % lo g( r i ! )- 

How should we modify this argument to deal with d- dimensional permu- 
tations? We fix a d- dimensional permutation X that is contained in A and 
consider a random ordering of all lines of the form A(i±, id, *)■ Given such 
an ordering, we use the chain rule to derive an upper bound on H(X). Each 
ordering yields a different bound. However, as in the one dimensional case, 
the key insight is that averaging over all possible orderings (in a class that 
we later define) gives us a simple bound on H(X). 

The overall structure of the argument remains the same. We consider a 
concrete line A(i\, id, *). Its contribution to the estimate of the entropy is 
log N where iV is the number of 1-entries that remain unshaded at the time 
(according to the chosen ordering) at which we compute the corresponding 
term in the chain rule for the entropy. However, now shade can fall from d 
different directions. The contribution of the line to the entropy will be the 
expected logarithm of the number of ones that remain unshaded after each 
of the d dimensions has cast its shade on it. 

The lines are ordered by a random lexicographic ordering. At the coarsest 
level lines are ordered according to their first coordinate i\. This ordering is 
chosen uniformly from § n . To understand how many l's remain unshaded 
in a given line, we first consider the shade along the first coordinate. If it 
initially has r 1-entries, then the number of unshaded 1-entries after this 
stage is uniformly distributed on [r]. We then recurse with the remaining 
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1-entries and proceed on the subcube of codimension 1 that is defined by the 
value of the first coordinate. It is not hard to see how the recursive expression 
for f(d, r) reflects this calculation. 

3.2 In detail 

Let A be a [n] d+1 -dimensional array of zeros and ones, and X is a random 
(i-permutation sampled uniformly from the set of d-permutations contained 
in A. Then H(X) = \og(Perd(A)) by theorem 2.2 and again we seek an 
upper bound on H(X). 

We think of X as an [n] d array each line of which contains each member 
of [n] exactly once. The proof does its accounting using lines of the form 
A(i\, ia, *), i.e., lines in which the (d+ l)-st coordinate varies. Such a line 
is specified by i = (ii,...,id). The random variable X\ is defined to be the 
value of X(ii, ...,id)- We think of the variables X\ as being revealed to us 
one by one. Thus, X ily4 _ iid must belong to 

Ri = Ri lt ...,i d = {j ■ A(i u ...,i d ,j) = 1} 

the set of 1-entries in this line. 

In the proof we scan these lines in a particular randomly chosen order. 
Let us ignore this issue for a moment and consider some fixed ordering of 
these lines. Initially, the number of 1-entries in this line is r\. As we proceed, 
some of these l's become unavailable to Xi, since choosing them would result 
in a conflict with the choice made in some previously revealed line. We say 
that these l's are in the shade of previously considered lines. This shade can 
come from any of the d possible directions. Thus we denote by Z x C i? ; the 
set of the indices of the 1-entries in R\ that are unavailable to X\ given the 
values of the preceding variables. We can express Z-^ = uf =1 Z^ where entries 
in Z\ are shaded from direction k. Namely, a member j of Ri belongs to Z\ 
if there is an already scanned line indexed by i' with X^ = j and where i 
and i' coincide on all coordinates except the k-th. Thus, given the values of 
the previously considered variables, there are at most 

iV; = \Ri \ Zi\ 

values that are available to Xi. 

We next turn to the random ordering of the lines. Now, however, we do 
not select a completely random ordering, but opt for a random lexicographic 
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ordering. Namely, we select d random permutations <Ti, a& G E> d - The line 
A{i\, id, *) precedes A(i[, ...,i' d ,*) if there is a A: G [n] such that o~k(ik) < 
a k{i'k) an d ij = i'j for all j < k. Thus a choice of the orderings induces a 
total order on the lines A(ii, i d , *). Denote this order by -<. That is, we 
write i — < j if i comes before j. We write i — <fc j if i — < j and i and j differ only 
in the A;-th coordinate. 

We think of X\ as being revealed to us according to this order. 

We turn to the definition of R\, Z\ and N\. Their definitions are affected 
by the chosen ordering of the lines. In addition, for reasons to be made 
clear later, we generalize the definition of N\. It is defined as the number 
of values available to X\ (given the preceding lines) from a given index set 
W C R{. In the discussion below, we fix X, a d-dimensional permutation 
that is contained in A. 

Definition 3.1. The index set of the 1-entries in the line A{i\, id, *) is 
denoted by 

Ri = Ri lr ..,i d = {j ■ A{i x , ...,i d ,j) = 1}, 
and its cardinality is = \Ri\. 

Let W C Ri with i = ...,id), and suppose that Xi e W. For a given 
ordering -<, let 

Z\[X, -<) = {j G Ri : X v = j for some V -< k i}. 

N i (W,X,^) = \W^U d k=1 Zl(X,^)\. 

Thus, Ni is a function of W C X and the ordering -<. Each variable 
X{ specifies a 1 entry of the line A(ii, id, *)■ The entry thus specified 
must conform to the values taken by the preceding variables. Namely, no 
line of X can contain more than a single 1 entry. We consider the number of 
values that the variable X\ can take, given the values that precede it. Fix an 
index tuple i = ...,i d ). The variable X± must specify an index with 
A(ii, i d +i) = 1, i.e., an element of R{. Consider some element j G R[. 
If Xi> = j, for some i' i and k < d then clearly Xi ^ j, or else the 
line X{ji, ik-i, *, ik+i, ■■■,id) contains more than a single j-entry. In other 
words, Xi cannot specify an element of Z\[X, -<) and is restricted to the set 
.Ri\Ufc =1 Z ; fc (X, -<). Therefore, there are at most N^Ri, X, -<) possible values 
that Xi can take given the variables that precede it in the order -<. 
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For a given order -<, we can use the chain rule to derive 
H(X) = J2H(X i \X r .)^i). 

i 

By theorem 2.2, 

: j -< i) = E Xj ^i [fr(Xi|Xj = xj : j -< i)] 

< Ejfj^i [log(^i(i2,, X, -<))] = E x [logCAT^^, X, -<))] . 

The last equality holds because N depends only on the lines of X that 
precede X\, and so taking the expectation over the rest of X doesn't change 
anything. 

As in the one dimensional case, the next step is to take the expectation 
of both sides of the above inequality over -<. 

H{X) <J2 E -< [MM^i, X, ■<))}} 

i 

= J2 E x [E-< (log(iV,(i2i,X } -0)]]. 

i 

The key to unraveling this expression is the insight that the random variable 
Ni has a simple distribution (as a function of -<), and moreover, that this 
distribution does not depend on X. 

Recall that in the one dimensional case, we obtained the distribution of 
Ni as follows. Initially, the number of ones in the i-th row was 7v Then the 
rows preceding the z-th row were revealed, and some of the ones in the i-th 
row became unavailable to X, because some other row had placed a one in 
their column. We defined Ni — \Ri \ Zi(a)\. The size of Zi(a) was shown to 
be uniformly distributed over {0, fj — 1}, and thus the distribution of N 
was shown to be uniform over {1, r^}. 

A similar argument works in the d dimensional case, but the distribution 
of Ni is no longer uniform. Recall that the function / is defined recursively 
by 

/(0,r) = log(r) 

Mr) = ^/(d-l,fc). 

r k=i 
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Claim 3.2. Let X be a d-permutation, i= (ii, ...,id) and let W C i?j be an 

index set such that Xi G W . Then E^ \log(Ni(W, X, -<))] depends only on d 
and r = \W\, and 

E^ [\og(N i (W,X,^))) = f(d,r). 

Proof. The proof proceeds by induction on d. 

First, note that if \W\ = r and d = 0, then Ni{W,X, -<) = \W\ = r by 
definition, and therefore 

E^ pog(JV,(Wi = log(r) = /(0, r). 

In order to proceed with the induction step, we must describe N\{W, X, -<) 
in terms of parameters of dimension d—1 instead of d. To this end we need 
the following definitions: 

• X' — X(ii, *,..., *). Note that X' is a (d — 1) -dimensional permutation. 

• W' = W \ -<). Note that | W| actually depends only on <Ti, the 
ordering of the first coordinate. 

• Let i' = = (i 2 , irf). 

• Given an ordering -<, let -<' be the ordering on the index tuples i' d _i) 
defined by the orderings 02, cr 3 , cr^- 

Note that for every i and -< we have JVi(W, X, -<) = N V (W, X', ^'). 

This equality follows directly from the definition of N. Now, 

E^ [log(Ar.(W, X, ■*<))] = E ai [E^, X, ^))]] 

= E CT1 [E K , [log(^(W, X', <))]] = [/(d - 1, \W'\)} 
The last step follows from the induction hypothesis. Consequently, 

E^ [\og(N i (W,X,^))] = J2P<\W'\ = k)f(d-l,k). 

k 

The only remaining question is to determine the distribution of \W'\ as a 
function of o~\. Note, however, that we have already answered this question in 
the one dimensional proof, namely, | W'\ is uniformly distributed on {1, r}. 
Indeed, W = |W \ Zl(X, -<)|, and Z\{X, -<) is the set of indices s such that: 
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• For some j G W , X(s, i 2 , id) = j (there are r — 1 such indices, one 
for each j e W). 

• The random ordering o\ places s before i\. 

In a random ordering, the position of %\ is uniformly distributed. Therefore 
\Z\(X, -<)| is uniformly distributed on {0, ...,r — 1}, and Pr(|W'| — k) — - 
for every 1 < k < r. 

Putting this together, we have shown that 

pog(JVi(W; X, -<))) — -^/(d— l,k) — f(d, r). 

k=l 

□ 

In conclusion, we have shown that 

H(X) <J2 E * [logiN^X, -<))]] 

i 

= J> x [/(d,n)] = X)/(d,n), 

i i 

where r- Y = r ilv . i(J is the number of ones in the vector A(ii, ^, *). There- 
fore, 

Per d {A) <]Je mn \ 

i 

4 The number of d- permutations — An upper 
bound 

As mentioned, the upper bound on the number of (/-dimensional permuta- 
tions is derived by applying theorem 1.6 to the all-ones array J . The main 
technical step is a derivation of an upper bound on the function f(d, r). 

Theorem 4.1. For every d there exist constants Cd and r d such that for all 
r > r d , 

log d (r) 

f(d,r)<hg(r)-d + c d ^^. 

r 

One possible choice that we adopt here is r d = e d for every d, c\ = 5, c 2 = 8, 
and Cd = d for d > 3 . 
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Proof. A straightforward induction on d yields the weaker bound f(d,r) < 
log(r) for all d, r. For d = there is equality and the general case follows 

since f(d,r) = ± Yll=i f( d ~ *) ^ 7ELi lo g( fc ) < log( r )- Tni s simple 
bound serves us to deal with the range of small r's (below r^-i). We turn to 
the main part of the proof 



rd-i 



fc=r d _i+l 



fc=l 



< 



log(r d _i) + ^ log(fe) - (d - 1) + c d _ 



fc=i 



loQfc) 



fe=l 

where £ = r^-i log(r^-i) = (d— l)e . It is easily verified that for r > > 3 
there holds log(r!) < rlog(r) — r + 21og(r). We can proceed with 



< - + log(r) H a H 



£ 

fc=i 



log^(fc) 
k 



We now bound the sum EI-=i log fc — by means of the integral J7 



log>) 



Note that the integrand is unimodal and its maximal value is 7 



Thus, 



c d ^j^ log d -\k) < c d _x ^ log d (r) | _ } 



r * — ' A; 

fc=i 



Putting this together, we have the inequality 



f(d, r) < log(r) — d + 



21og(r) + £ + c,_ 1 ( 7 + ^ 



Therefore it is sufficient to choose q such that for every r > e a 



2 log(r) + £ + c d _! ( 7 + ) < c d log d (r) 
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i.e., 



log 



d-l, 



+ 



log d (r) 



Cd-l 



1 



log>) 



1 



The left hand side of the above inequality is clearly a decreasing function of 
r. Therefore it is sufficient to verify the inequality for r = e d . Plugging this 
and the values of the constants £ and 7 into the left hand side of the above 
inequality, we get 



d d - 1 
< 

Thus, we may take 
c d = 



(d-l)e 
1 

e d-i 
1 



d-l 



+ C d -l 



(d 



\d-l 



D d-1 



d d 



Cd-l 



Cd-l 



i 4 



d d 1 \d) 



ey 



od-l 



d 



ey 
d) 



Calculating Cd using this recursion and the fact that cq = 0, we get that 

ci = 2 + . 
induction, 



c\ = 2 + e < 5, C2 < 8, and q < d for 3 < d < 10. Proceeding by 



c d 



1 



1 



od-l 



(d- 1) 3 (1.1) 



+ d 



< 



[i.iY(d 



d\ 



+ 2d(-) < 



;i.i) d (d 



2 /eN 



2c/ 2 (l.l)<*d 3 



d! d\ 

In the inequality before the last one, we used the fact that for d > 10 , 

\d) — d\- 

□ 

For the [n] d+l all ones array J, r^.,,^ = n for every tuple (ii, id), and 
so for large enough n we have the bound 



Per d (J)< ' e 



f(d,n) 



(e^) n < (exp 



Hi— M 



log(n) -d + c d 



\og\n) 



n 



For a constant d, letting n go to infinity, Cd log ^ 
number of d permutations is at most 



oil) and therefore the 



16 



References 

[1] L. M. Bregman, Certain properties of nonnegative matrices and their 
permanents, Dokl. Akad. Nauk SSSR 211 (1973), 27-30. MR MR0327788 
(48 #6130) 

[2] T. M. Cover and J. A. Thomas, Elements of Information Theory, 
Wiley, New York, 1991. 

[3] G.P. Egorichev, Proof of the van der Waerden conjecture for perma- 
nents, Siberian Math. J. 22 (1981), 854-859. 

[4] D.I. Falikman, A proof of the van der Waerden conjecture regarding 
the permanent of a doubly stochastic matrix, Math. Notes Acad. Sci. 
USSR 29 (1981), 475-479. 

[5] L. Gurvits , Van der Waerden /SchrijverValiant like conjectures and 
stable (aka hyperbolic) homogeneous polynomials: one theorem for all. 
With a corrigendum, Electron. J. Combin. 15 (2008), R66 (26 pp). 

[6] M. Jerrum, A. Sinclair, and E. Vigoda, A polynomial-time ap- 
proximation algorithm for the permanent of a matrix with nonnegative 
entries, J. ACM, 671-697. 

[7] M. Kochol, Relatively narrow Latin parallelepipeds that cannot be ex- 
tended to a Latin cube, Ars Combin. 40 (1995), 247260. 

[8] Laurent, M. and Schruver, A., On Leonid Gurvitss proof for per- 
manents, The American Mathematical Monthly 10 (2010), 903-911. 

[9] H. Minc, Permanents, Encyclopedia of Mathematics and Its Applica- 
tions Vol. 6, Addison- Wesley, Reading, Mass, 1978 

[10] Jaikumar Radhakrishnan, An entropy proof of Bregman's theorem, 
J. Combinatorial Theory Ser. A 77 (1997), no. 1, 80-83. MR MR1426744 
(97m: 15006) 

[11] Schruver, A., A short proof of Mine's conjecture, J. Comb. Theory 
Ser. A 25 (1978), 80-83. 

[12] L. G. Valiant, The complexity of computing the permanent, Theoret. 
Comput. Sci. 8 (1979), 189-201 



17 



J. H. van Lint and R. M. Wilson, A Course in Combinatorics, 
Cambridge U.P., 1992. 

B.D. McKay and I.M. Wanless, A census of small latin hypercubes, 
SIAM J. Discrete Math. 22 (2008), pp. 719736. 



18 



